STABILITY OF THE SLOW MANIFOLD 
IN THE PRIMITIVE EQUATIONS 
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Abstract. We show that, under reasonably mild hypotheses, the solution of 
the forced-dissipative rotating primitive equations of the ocean loses most of 
its fast, inertia— gravity, component in the small Rossby number limit as t — > 
oo. At leading order, the solution approaches what is known as "geostrophic 
balance" even under ageostrophic, slowly time-dependent forcing. Higher- 
order results can be obtained if one further assumes that the forcing is time- 
independent and sufficiently smooth. If the forcing lies in some Gevrey space, 
the solution will be exponentially close to a finite-dimensional "slow manifold" 
after some time. 



1. Introduction 

One of the most basic models in geophysical fluid dynamics is the primitive 
equations, understood here to be the hydrostatic approximation to the rotating 
compressible Navier-Stokes equations, which is believed to describe the large-scale 
dynamics of the atmosphere and the ocean to a very good accuracy. An important 
feature of such large-scale dynamics is that it largely consists of slow motions in 
which the pressure gradient is nearly balanced by the Coriolis force, a state known as 
geostrophic balance. Various physical explanations have been given, some supported 
by numerical simulations, to describe how this comes about, but to our knowledge 
no rigorous mathematical proof has been proposed. (For a review of the geophysical 
background, see, e.g., [7j.) One aim of this article is to prove that, in the limit 
of strong rotation and stratification, the solution of the primitive equations will 
approach geostrophic balance as t — > oo, in the sense that the ageostrophic energy 
will be of the order of the Rossby number. 

As illustrated by the simple one-dimensional model (|4.3j) . here the basic mech- 
anism for balance is the viscous damping of rapid oscillations, leaving the slow 
dynamics mostly unchanged. Separation of timescale, characterised by a small pa- 
rameter £, is therefore crucial for our result; this is obtained by considering the 
limit of strong rotation and stratification, or in other words, small Rossby number 
with Burger number of order one. We note that there are other physical mecha- 
nisms through which a balanced state may be reached. Working in an unbounded 
domain, an important example is the radiation of inertia-gravity waves to infinity 
in what is known as the classical geostrophic adjustment problem (see §7.3] 
and further developments in [27j). 
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Attempts to extend geostrophic balance to higher orders, and the closely related 
problem of eliminating rapid oscillations in numerical solutions (e.g., [3l l20lfT6l l34j). 
led naturally to the concept of slow manifold j!9j . which has since become impor- 
tant in the study of rotating fluids (and more generally of systems with multiple 
timescales). We refer the reader to [21] for a thorough review, but for our pur- 
poses here, a slow manifold means a manifold in phase space on which the normal 
velocity is small; if the normal velocity is zero, we have an exact slow manifold. 
In the geophysical literature, there have been many papers proposing various for- 
mal asymptotic methods to construct slow manifolds (e.g., [351 13Z]). A number of 
numerical studies closely related to the stability of slow manifolds have also been 
done (e.g., P21E5]). 

It was realised early on [19l [36] that in general no exact slow manifold exists and 
any construction is generally asymptotic in nature. For finite-dimensional systems, 
this can often be proved using considerations of exponential asymptotics (see, e.g., 
|15|). More recently, it has been shown explicitly [33] in an infinite-dimensional 
rotating fluid model that exponentially weak fast oscillations are generated spon- 
taneously by vortical motion, implying that slow manifolds could at best be ex- 
ponentially accurate (meaning the normal velocity on it be exponentially small). 
Theorem [5] shows, given the hypotheses, that an exponential accuracy can indeed 
be achieved for the primitive equations, albeit with a weaker dependence on e. 

From a more mathematical perspective, our exponentially slow manifold (see 
Lemma [5]) , which is also presented in [31j in a slightly different form, is obtained 
using a technique adapted from that first proposed in 22J. It involves truncating 
the PDE to a finite-dimensional system whose size depends on e and applying 
a classical estimate from perturbation theory to the finite system. By carefully 
balancing the truncation size and the estimates on the finite system, one obtains 
a finite-dimensional exponentially accurate slow manifold. This estimate is local 
in time and only requires that the (instantaneous) variables and the forcing be in 
some Sobolev space H s ; it (although not the long-time asymptotic result below) 
can thus be obtained for the inviscid equations as well. If our solution is also Gevrey 
(which is true for the primitive equations given Gevrey forcing), the ignored high 
modes are exponentially small, so the "total error" (i.e. normal velocity on the slow 
manifold) is also exponentially small. 

Gevrey regularity of the solution is therefore crucial in obtaining exponential 
estimates. As with the Navier-Stokes equations [9], in the absence of boundaries 
and with Gevrey forcing, one can prove that the strong solution of the primitive 
equations also has Gevrey regularity [25] . For the present article, we need uniform 
bounds on the norms, which have been proved recently [23] following the global 
regularity results of [6] [T3J [14]. Since our result also assumes strong rotation, 
however, one could have used an earlier work [2] which proved global regularity 
under a sufficiently strong rotation and then used [25] to obtain Gevrey regularity. 

While our earlier paper [31j is concerned with a finite-time estimate on pointwise 
accuracy ( "predictability" ) , in this article our aim is to obtain long-time asymptotic 
estimates (on "balance"). In this regard, the main problem for both the leading- 
order (Theorem!]} and higher-order (Theorem[2]) estimates are the same: to bound 
the energy transfer, through the nonlinear term, from the slow to fast modes at 
the same order as the fast modes themselves. For this, one needs to handle not 
only exact fast-fast-slow resonances, whose absence has long been known in the 
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geophysical literature (cf. e.g., [U [8] \17\ [35] for discussions of related models), but 
also near resonances. A key part in our approach is an estimate involving near 
resonances in the primitive equations (cf. Lemma [1} . Another method based on 
algebraic geometry to handle related near resonances can be found in pQ. 

Taken together with [31] . the results here may be regarded as an extension of 
the single-frequency exponential estimates obtained in }22j to the ocean primitive 
equations, which have an infinite number of frequencies. Alternately, one may view 
Theorem [5] as an extension to exponential order of the leading-order results of [2 
for a closely related model. Finally, our results here put a strong constraint on the 
nature of the global attractor [12] in the strong rotation limit: the attractor will 
have to lie within an exponentially thin neighbourhood of the slow manifold. 

The rest of this article is arranged as follows. We begin in the next section by 
describing the ocean primitive equations (henceforth OPE) and recalling the known 
regularity results. In Section[3l we write the OPE in terms of fast-slow variables and 
in Fourier modes, followed by computing explicitly the operator corresponding to 
the nonlinear terms and describing its properties. In Section 0] we state and prove 
our leading-order estimate, that the solution of the OPE will be close to geostrophic 
balance as t — > oo. In the last section, we state and prove our exponential-order 
estimate. 

2. The Primitive Equations 

We start by recalling the basic settings of the ocean primitive equations [TS] . 
and then recast the system in a form suitable for our aim in this article. 

2.1. Setup. We consider the primitive equations for the ocean, scaled as in [24] 

d t v + - [v 1 - + W 2 p] + u-Vv = pAv + f v> 
1 , 

(2.1) d t p- -u +u-Vp = pAp + f p , 

V • u = V-v + d z u s = 0, 
p = -d z p. 

Here u = (it 1 , u 2 , u 3 ) and v — (u l ,u 2 ,0) are the three- and two-dimensional fluid 
velocity, with v := (— it 2 , v}, 0). The variable p can be interpreted in two ways: 
One can take it to be the departure from a stably-stratified profile (with the usual 
Boussinesq approximation), with the full density of the fluid given by 

(2.2) pw,(x, y, z, t) = po - e~ x zpi + p(x, y, z, t), 

for some positive constants po and p\. Alternately, one can think of it to be, 
e.g., salinity or temperature that contributes linearly to the density. The pressure 
p is determined by the hydrostatic relation d z p = —p and the incompressibility 
condition V-tt = 0, and is not (directly) a function of p. We write V := (d x ,d y , d z ), 
V 2 := (d x ,d y , 0), A := 9 2 +£? 2 + i9 2 and A 2 := 9 2 +<9 2 . The parameter e is related to 
the Rossby and Froude numbers; in this paper we shall be concerned with the limit 
e — > 0. In general the viscosity coefficients for v and p are different; we have set 
them both to p for clarity of presentation (the general case does not introduce any 
more essential difficulty). The variables (y, p) evidently depend on the parameters 
e and p as well as on (x,t), but we shall not write this dependence explicitly. 
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We work in three spatial dimensions, x := (x,y,z) = (x 1 ,x 2 ,x 3 ) £ [0,Li] x 
[0,1(2] x [—L3/2, £3/2]=: M, with periodic boundary conditions assumed; we write 
|M| := L1L2L3. Moreover, following the practice in numerical simulations of strat- 
ified turbulence (see, e.g., [4]), we impose the following symmetry on the dependent 
variables: 

v(x,y,-z) = v(x,y,z), p(x,y,-z) = p(x,y, z), 

(2.3) 

u 3 (x, y, -z) = -u 3 (x, y, z), p(x, y, -z) = -p(x, y, z); 

we say that v and p are even in z, while u 3 and p are odd in z. For this symmetry to 
persist, f v must be even and / odd in z. Since u 3 and p are also periodic in z, we 
have u 3 {x, y, -L 3 /2) = u 3 (x, y, L 3 /2) = and p(x, y, -L 3 /2) = p(x, y, L 3 /2) = 0; 
similarly, dzU 1 = 0, d z u 2 — and d z p = on z — 0, ±L 3 /2 if they are sufficiently 
smooth (as will be assumed below). One may consider the symmetry conditions 
(|2.3|) as a way to impose the boundary conditions u 3 = 0, p — 0, d z u x = 0, <9jW 2 = 
and d z p = on z = and z — L 3 /2 in the effective domain [0, Li] x [0, L2] x [0, £3/2]. 
All variables and the forcing are assumed to have zero mean in M; the symmetry 
conditions above ensure that this also holds for their products that appear below. 
It can be verified that the symmetry (|2.3[) is preserved by the OPE (|2 . lj) : that is, 
if it holds at t = 0, it continues to hold for t > 0. 

2.2. Determining the pressure and vertical velocity. Since u 3 = at z = 0, 
we can use (|2.1b 'l to write 

(2.4) u 3 (x,y,z) = - [ V-v(x,y,z')dz'. 



Similarly, the pressure p can be written in terms of the density p as follows (cf. 
|28j). Let p(x,y,z) = (p(x,y)) + 5p(x,y, z) where (•) denotes z-average and where 



(2.5) 5p{x,y,z) = - p(x,y,z')dz' 

J zo 

with zq(x, y) chosen such that (Sp) = 0; this is most conveniently done using Fourier 
series (see below). Using the fact that 

(2.6) / V-vdz = - d z u 3 dz = u 3 (-,-L 3 /2) - u 3 (-,L 3 /2) = 0, 
J-L3/2 V-L3/2 

and taking 2d divergence of the momentum equation (|2.1a ). we find 

(2.7) -[V ■ (v^ + A 2 (p)] + V-(u-Vv)=pAV-(v)+V-(f v ). 

e 

Here we have used the fact that z-integration commutes with horizontal differential 
operators. We can now solve for the average pressure (p), 

(2.8) (p) = A^ 1 [-V • (v x ) + e(-V • (u ■ Vv) + /iAV ■ (v) + V • (/„))] 

where A^ 1 is uniquely defined to have zero xy-average. With this, the momentum 
equation now reads 

(2 9) dtV + e ^ V± ~ VA ^ ly ' ( v± ^ + ^ 2& p\ + u • - VA2 X V • (u ■ V«) 

= pA(v - VA 2 - J V • («)) +/. - VA2 X V • (f v ). 
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2.3. Canonical form and regularity results. Besides the usual L P (M) and 
iP(M), with p £ [l,oo] and s > 0, we shall also need the Gevrey space G" T (M), 
defined as follows. For a > 0, we say that u £ G" T (M) if 

(2.10) |e ff( - A)1/2 u| L2 =: M G . < oo. 

Let us denote our state variable W = (v, p) T . We write W £ L P (M) if v £ L P (M) 2 , 
p G L P (M), (v, p) has zero average over M and (v, p) satisfies the symmetry (|2.3[) . in 
the distribution sense as appropriate; analogous notations are used for W £ H S (M) 
and W £ G CT (M), and for the forcing / (which has to preserve the symmetries of 
W). 

With u 3 given by (|2^|) and 5p by (|2~5j) . we can write the OPE f2TTb) and (|23|) 
in the compact form 

(2.11) a t PF+iLVK + S(VK,VF)+AVK = /. 
The operators L, B and A arc defined by 

LVF = (v 1 - - VAJ X V ■ (t^) + V 2 5p, -w 3 ) T 

(2.12) B(W, #) = (it • V* - VA^V • (u ■ Vu},« ■ Vp) T 

AW = -(a*A(« - VA^V • («)), pAp) T , 
and the force / is given by 

(2.13) f=(f v -VA- 1 V-(f v )J p ) T . 

The following properties are known (see, e.g., 25J). The operator L is antisymmet- 
ric: for any W £ L 2 (M) 

(2.14) (LW,W) L 2=0: 

B conserves energy: for any W £ H 1 (JA) and W £ H 1 (3vV), 

(2.15) (W,B(W,W)) L 2 = 0; 

and A is coercive: for any W £ H 2 (JA), 

(2.16) (AW, W) L 2 = (j, \VW\ 2 L2 . 

We shall need the following regularity results for the OPE (here K s and M a are 
continuous increasing functions of their arguments): 

Theorem 0. Let W £ H 1 and f £ L°°(IR + ; L 2 ). Then for allt>0 there exists a 
solution W(t) £ H 1 of (LTTDl with W(0) = W and 

(2.17) \W(t)\ m <Jfo(|Wb| ff i,||/||o) 

where, here and henceforth, \\f\\ s ■= esssup t>0 \f(t)\ H s for s > 0. Moreover, there 
exists a time Ji(|Wo|iji , ||/||o) such that for t > T\, 

(2.18) |W(t)| Hl < Xidl/y. 

Similarly, if f £ H s ~ x ), there exists a time T s (\Wq\hi, such that 

(2.19) mt)\ B .<K.(\\f\\ a _ 1 ) 

for t > T s . Finally, fixing a > 0, if also V/ £ L°°(R+;G" T ) ; there exists a time 
T CT (|Wo| H1 , |V/| G such that, for t > T a 

(2.20) \V 2 W(t)\ G „ <M a (\Vf\ G «). 
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The proof of (pTTTjl — ([2TT8)) can be found in [12]; the higher-order results ([27T9]) can 
be found in [23]. Both these works followed [6] and [13]. The result ()2.20p follows 
from [35] and using (|2.19|) for s = 2. 

Since we are concerned with the limit of small e, however, one might also be 
able to obtain (|2.17[) and (]2 . 19[) following the method used in [2] for the Boussinesq 
(non-hydrostatic) model. One could then proceed to obtain (|2.20p as above. 

3. Normal Modes 

In this section, we decompose the solution W into its slow and fast components, 
expand them in Fourier modes, and state a lemma that will be used in sections [4] 
and E] below. 

3.1. Fast and slow variables. The Ertel potential vorticity 

(3.1) q E =V ± -v-d zP + e [(d z v) ■ V ± p - d zP (V ± -v)] , 

where V := (—d y , d x , 0), plays a central role in geophysical fluid dynamics since it 
is a material invariant in the absence of forcing and viscosity. In this paper, however, 
it is easier to work with the linearised potential vorticity (henceforth simply called 
potential vorticity) 

(3.2) q:=V^-v-d z p, 

From (|2.1|) . its evolution equation is 

(3.3) d t q + V ± -{u- Vv) -d z (u- Vp) = pAq + f q 

where f q := V -f v ~ d z f p . Let ■0° :— A^ 1 ^, uniquely defined by requiring that ip° 
has zero integral over M, and let 



(3.4) : ~ \p° J '~ \-d z ip° 

We note a mild abuse of notation on v° and V -1 : W° = (—d y ip°, d x ip°, —d z ip°) T . 

A little computation shows that W° lies in the kernel of the antisymmetric 
operator L, that is, LW° = 0. Conversely, if LW = 0, then W = (V x *, -~d z ^) T 
for some Since u 3 = 0, we have V-u = 0, so v = V ±, 3 / + V for some ^(x, y, z) 
and V(z). Now 

= v ± - V 2 A^ 1 V-(v- L ) + V 2 Sp 

(3 5) 

= -V 2 * + F x + V 2 A2 1 A 2 (#) + V2<5p. 

Since all other terms are horizontal gradients and V does not depend on (x,y), we 
must have V = 0. Writing "i>(x,y, z) — ^>(x,y, z) + (^(x,y)) where ^>(x,y,z) has 
zero 2-average, the terms that do not depend on z cancel and we are left with 

(3.6) -V 2 * + V 2 5p = 0. 

So Sp(x, y, z) — ^f(x, y, z) + $(z); but since (Sp) = 0, $ = and thus p = —d z ^ by 
(|2.5|) . Therefore the null space of L is completely characterised by (|3.4I) . 

(3.7) kerL = {W° : W° = (V x ^°, -d z iJj ) T }. 

With tp° = A _1 (V ± -v — d z p) as above, this also defines a projection W i— > W°. 
We call W° our slow variable. 
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Letting B° be the projection of B to kerL, 

/ V ± A- 1 [V- 1 • (u ■ Vi>) - dJu ■ Vp)l 

we find that W° satisfies 

(3.9) d t W° + B°(W,W) + AW° = f 

where f° = (V^A -1 /^, —d z A~ 1 f q ) T is the slow forcing. 
Now let 

(3.10) W '=(£y.=w-w°=(j t zpy 

It will be seen below in Fourier representation that W £ is a linear combination of 
eigenfunctions of L with imaginary eigenvalues whose moduli are bounded from 
below; we thus call W £ our fast variable. Since V-i>° = 0, the vertical velocity u 3 
is a purely fast variable. In analogy with (|3.9p . we have 

(3.11) d t W £ + -LW £ +B E (W, W) + AW £ = f 

e 

where B £ {W, W) := B{W, W) - B°(W, W) and f £ := f - f°. 

The fast variable has no potential vorticity, as can be seen by computing V • 
v £ — d z p £ = q — V^-V^V — d zz ip° = 0. Since the slow variable is completely 
determined by the potential vorticity, this implies that the fast and slow variables 
are orthogonal in L 2 (M), 

(W°, W £ ) L * = (v°,v £ ) L 2 + (p°, p £ ) L , 

' ' ' = (V V, v £ ) L2 - (d z ^°, P £ ) L 2 = (v> , -v x • v £ + 8 zP £ ) L 2 = 0. 

Of central interest in this paper is the "fast energy" 

(3.13) ||W £ |l 2 = i(K|l 2 + |p £ H 2 ). 

Its time derivative can be computed as follows. Using (|3.12|) , we have after inte- 
grating by parts 

(3.14) (W £ ,d t W) L * = (W £ ,d t W°) L 2 + (W £ ,d t W £ ) L 2 = ±±-\W £ \ 2 L2 . 
Now P7T5|) implies that 

(3.15) {W S ,B(W,W)) L 2 = (W^iW^W + W £ )) L 2 = {W £ ,B{W 1 W°)) L 2. 
Putting these together with (|2. 14|) and (|2.16|) , we find 

(3.16) \^\W £ \ 2 L 2 + p\VW £ \ 2 L 2 = -{W S ,B(W,W°)) L 2 + (W £ J £ ) L 2. 

3.2. Fourier expansion. Thanks to the regularity results in Theorem 0, our so- 
lution W(t) is smooth and we can thus expand it in Fourier series, 

(3.17) v(x,t)=j: k v k (t)e ikx and p(x, t) = J2 k Pk (t) c ikx . 

Here k = (ki,k 2 ,k 3 ) € Z L where Z L = R 3 /M = {(2wh/Li, 2nl 2 /L 2 , 27rZ 3 /L 3 ) : 
(h, h, h) S ^ 3 }; any wavevector k is henceforth understood to live in Z^. We also 
denote k' :— (ki,k 2) 0) and write k' A j := k\j 2 — k 2 j\. Since our variables have 
zero average over M, v k = when k = 0; moreover, since p is odd in z, p k = 
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whenever k 3 = 0. Thus Wk '■= (vk,pk) = when k = 0, which allows us to write 
the H s norm simply as 

(3.18) \W\ 2 H ,=Ek\k\ 2s \Wk\ 2 

and (see (l2~10|) for the definition of G a ) 

(3.19) \W\ 2 G ,=J2k^ lkl \Wk\ 2 . 

The antisymmetric operator L is diagonal in Fourier space, meaning that L k i = 
when k jtz I; we shall thus write L k := L k k- When k 3 ^= 0, we have 

/ -1 -fci/fc 3 N 

(3.20) L fe = 1 -k 2 /k 3 

\h/k 3 k 2 /k 3 

For k' 0, its eigenvalues are 0^ = and iwj^ = ±i|fc|/fc3, where |fc| := (k\ + k. 
fcf) 1 / 2 , with eigenvectors 

/ fc 2 \ i /-k 2 k 3 ±ik 1 \k\ 

(3 . 21) X 2 = m (-£ j a„d X t = ^ Mj± U, W 

When fc' = 0, we have u k = and iw^ = ±i as eigenvalues with eigenvectors 

(3.22) 1^=1 I and 



2 





For k fixed, these eigenvectors are orthonormal under the inner product • in C 3 . 

When k 3 = 0, the fact that pk = and k ■ Vk — implies that the space is one- 
dimensional for each k (in fact, it is known that the vertically-averaged dynamics 
is that of the rotating 2d Navier-Stokes equations). Since projecting to the k 3 = 
subspace is equivalent to taking vertical average, we compute 

(3.23) (LW) = ((it 1 -) - V a Aa 1 V-<ir L },0) T 

where we have used (it 3 ) — (since it 3 is odd) and (dp) = (by definition). 
Reasoning as in (|3.5p - (|3.6p above, we find that (LW) = 0, that is, the vertically- 
averaged (k 3 — 0) component is completely slow. In this case we can thus write 



2 



1 / k 

(3.24) L0° k = and X k = Jj^ 

which can be included in the generic case k' ^ in computations. Since the k 3 = 
component is completely slow, (W e ) = 0, there is no need to fix 

J k 



We note that since k 3 ^ 0, \u> k \ > 1, viz. 



,±|2 _ m £ / ^1 + fej + ^3 



(3.25) inf = inf i ' 3 , 1=1. 



In what follows, it is convenient to use {X°,X±} as basis. 
We can now write 

W°{x,t) :=£ fc w° k {t)X°e ikx 
(3.26) . . , ., 

W e (x,t) := Efe ^IW^e-'^/V^, 
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where s € {— 1, +1}, which we write as {— , +} when it appears as a label. The 
Fourier coefficients w k and w k are complex numbers that depend on t only, with 



q = and w^ ki k ^ Q . = 0. With a 6 { — 1, 0, +1}, they can be computed using 

(3.27) w% (t) = ^ j W(x , t) ■ XI ds. 
The following relations hold: 

(3.28) \W*\\,=Y. k \<? and \W%, = £* K| 2 - 
In addition, the fact that (v°, p°) is real implies 

(3.29) w°_ k = -wl and wf klM _ fcs) = w^MM) 

where overbars denote complex conjugation. Similarly, since (v E ,p 6 ) is real, 

(3.30) w ± k = w^ and w % u k 2 ,-k 3 ) = ~ w %iMM) 
when k! ^ and, when fe' = 0, 

(3.31) wf 00 _ k3) =wT Qk3y 

We shall see below that, the linear oscillations having been factored out, the variable 
w k is slow at leading order. Similarly to W, we write the forcing / as 

f(^):=E;/*(i)^, 

where, unlike in (|3.26|) . there is no factor of e -IWfct / £ in the definition of f £ . As 
noted above, / must satisfy the same symmetries as W , so the above properties of 



w k also hold for f k ; we note in particular that f k = when /c 3 = 0. 
For later convenience, we define the operator <9 t * by 

(3.33) d* t W := e- tL ' e d t e tL/s W. 
From (|2.11|) . we find 

(3.34) dfW + B(W,W) + AW = f, 

which is dfW with the large antisymmetric term removed. 

Now the nonlinear term on the rhs of (|3.16|) can be written as 

{W £ , B(W° + W £ , W°)) L 2 = (W £ , B(W°, W°)) L 2 + {W e ,B(W £ ,W°)) L 2 

(3 ' 35) = {W £ ,B(W°,W°)) L 2 - (W°,B(W £ ,W £ )) L 2, 

where the identity (W°, B(W £ , W £ )) L 2 = -(W £ , B(W £ , W°)) L 2 had been obtained 
from 

First, let 

(W £ , B(W°, W°)) L 2 = |M| w Q jW ° k ^i(X° ■ k')(X° k ■ XI) 6 j+k - t e^ £ 

(3.36) s jkl 

w j w k w l B jkl e ' 

jkl 

where Sj+ k -i = 1 when j + k = I and otherwise, and where 

(3.37) Bf k i := i |M| 6 j+k . t (X9 ■ k')(X° k ■ Xf). 
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It is easy to verify from ([37371) that B ^ = when \j'\ \k'\ h = 0, so we consider 
the other cases. For the first factor, we have 

(3.38) X^-k' = ^f. 

\3\ 

For the second factor, we have 

x o . X s = k 2 -isk 1 , = 

k 1 y/2\k\ 

V^|l!||i||('| 
From these, we have the bound 

(3.40) iBsi<«aa. 

Next, we consider 

(W°,B{W e ,W £ )) L 2 



(3.41) 



where 

(3.42) := i |M| «5 i+fe _i (VXr ■ fe)(X| • X°) 

and where the operator V, which produces an incompressible velocity vector out of 
Xj, is defined by 

VXJ = X} when j 3 \f\ = 0, and 

(3.43) wTrr 1 /"-J2j3+irji|j|^ 

V2b'lb"l 

Thus, we have VXT- • fc = when J3 = 0, 

(3.44) VX; ■ fe = (fci - irfc 2 ) /v/2 
when j = 0, and 

(3 45) VJT k - Mj ' A k ' ] + irb1(j/ ' fc,) ~ ir|j ' /|2|il h/j3 

V2\j\\j'\ 

in the generic case js|j | 7^ 0. In all cases, we have the bound 
(3.46) |VXJ ■ fe| < |M| (V2|fe'| + \j'\ |fc 3 |/|j 3 |). 



V ^ = 77^77777^7 I ^jjj^' 1 J when j3|/l * °' 
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Next, X s k ■ A° = when fc 3 = or k' = l' = 0, and 



Xt . X? = h t ish when I 1 £ and k' = 0, 

k 1 y/2\l\ 

\k'\ 

(3.47) X% ■ X? = sgal 3 -j=-^ when «' = and feV 0, 

! = Hk'-i')k 3+ is(k'An\ki + \k'i% when|fe , ||zWo . 

These give us the bound 

(3.48) \Xi-X?\<y/5/2 

in all cases and, together with Q3.46p , when % ^ 0, 

(3.49) l^l^v^lMKlfc'l + lj'Mfcal/ljal). 
When J3/C3 = or I = 0, we have Bjf^ = 0. 

3.3. Fast— Fast— Slow Resonances. We first write (|3.4ip as 

(3.50) (W°,B(W S , W*)) L * = \ f; + B&) e -^J+-J)tA. 

It has long been known in the geophysical community that many rotating fluid 
models "have no fast-fast-slow resonances" (see, e.g., [35] for the shallow- water 
equations and [4] for the Boussinesq equations). In our notation, the absence of 
exact fast-fast-slow resonances means that B^f^ + -Bjyj = whenever u>j +cuf. = 0; 
the significance of this will be apparent below [see the development following (|4.16p ] . 
For our purpose, however, we also need to consider near resonances, i.e. those cases 
when |<x^-|-ci;|.| is small but nonzero. The following "no- resonance" lemma contains 
the estimate we need: 

Lemma 1. For any j, k, I G with I 7^ 0, 

(3.51) + < c nr |M| (ffl + + N) K 
where c nr is an absolute constant. 

We note that = Sj^o = when I = by (|3~T1) . so this case is trivial. We 

defer the proof to Appendix [A] 

4. Leading-Order Estimates 

In this section, we discuss the leading-order case of our general problem. This 
is done separately due to its geophysical interest and since it requires qualitatively 
weaker hypotheses. As before, W(t) = W°(t) + W E (t) is the solution of the OPE 
(|2.11|) with initial conditions W(0) = Wo, and K s (-) is a continuous and increasing 
function of its argument. 

Theorem 1. Suppose that the initial data Wq £ H (M.) and that the forcing f € 
L 0O (R + ;iJ 2 )nM /1 ' 0O (M + ;i 2 ), with 

(4.1) ||/|| g := esssu Pt>0 (\f(t)\ H2 + \d t f(t)\ L2 ). 
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Then there exist T g = T g (| Wol^i , ||/|| g , e) and K g = K g (\\f\\ g ), such that for t > T g , 
(4.2) \W%t)\ L2 <V^K g (\\f\\ g ). 

In geophysical parlance, our result states that, for given initial data and forcing, 
the solution of the OPE will become geostrophically balanced (in the sense that 
the ageostrophic component W £ is of order y/e) after some time. We note that 
the forcing may be time-dependent (although ||/|| g cannot depend on e) and need 
not be geostrophic; this will not be the case when we consider higher-order balance 
later. Also, in contrast to the higher-order result in the next section, no restriction 
on e is necessary in this case. 

The linear mechanism of this "geostrophic decay" may be appreciated by mod- 
elling (|3.1ip . without the nonlinear term, by the following ODE 

dx i 
4.3 +~ x + fl x = f 

dt e 

where fx > is a constant and / = f(t) is given independently of e. The skew- 
herniitian term ix/e causes oscillations of x whose frequency grows as e — > 0. In 
this limit, the forcing becomes less effective since / varies slowly by hypothesis 
while the damping remains unchanged, so x will eventually decay to the order of 
the "net forcing" \fef . More concretely, let z(t) = e lt ^ e x{t) and write (|4.3p as 

(4.4) ^ (e^ 2 z) + |e^ 2 z = e^-^f, 
from which it follows that 

(4.5) ^(e^ 2 |z| 2 ) + ^/ 2 \z\ 2 = 2e^ 2 Re (e^zf). 
Integrating, we find 

3 M*/2| 



(4.6) 



\z{t)\ - |z(0)| 2 +fx [ e^ 2 |z(r)| 2 dr = 2 f e^ r/2 Re( e - ir/£ z/) dr 
Jo Jo 
ft 

= 2e[e^ T/2 Re(ie- it/e z/)]J ) -2£ / Re[ie- iT / £ 9 r (e^ T / 2 z/)] dr, 
"' Jo 



where the second equality is obtained by integration by parts. Since dtf is bounded 
independently of e, the integral can be bounded using (|4.4p and the integral on the 
left-hand side. This leaves us with 

(4.7) |z(t)| 2 < e"^ 2 Cl (|/|) |z(0)| 2 + f (1 - e-^ 2 ) K(\f\, |9 t /|, M ). 

Most of the work in the proof below is devoted to handling the nonlinear term, 
where particular properties of the OPE come into play. A PDE application of this 
principle can be found in [29j . 

4.1. Proof of Theorem [U In this proof, we omit the subscript in the inner prod- 
uct (•; Ol 2 when the meaning is unambiguous; similarly, | • | = | • L 2 . We start by 
writing (|3.16[) as 

^-\W e \ 2 + 2fx\VW e \ 2 
dt 

(4-8) = _ 2 (V7 £ , B(W°, W )) - 2{W E , B{W e , W a )) + 2(W e , f £ ) 

= -2{W e , B(W°, W )) + 2(W°, B(W £ , W e )) + 2{W £ , f). 



STABILITY OF THE SLOW MANIFOLD 



13 



Using the Poincare inequality, \W £ \ 2 < c p |VVU £ | 2 , and multiplying the left-hand 
side by 2e vt where v := /ic p , we have 

(4.9) -^(e^W^I 2 ) + ne ut \VW e \ 2 < e vt (^-\W e \ 2 + ^\VW £ \ 2 + /j|V1U £ | 2 ) . 
With this, P~5]) becomes 

^_( e "*|^| 2 ) +fxe ut \VW e \ 2 

(4.10) di v 1 1 1 P 1 1 

< 2e"* (VU £ , / £ ) - 2 e"* (J¥ £ , £(VF°, IU )) + 2e vt (W°, B(W £ , W £ )). 
We now integrate this inequality from to t. On the left-hand side we have 

(4.11) 



I ( T" (e^ I ^ 1 2 ) + Me" T | V 1U £ | 2 } dr 
.11) 7oldT J 



= e l "|^ £ (t)| 2 -|VK e (0)| 2 + M / e VT \VW e \ 2 dr. 

Jo 

Using the expansion (|3 . 26[1 of 1U £ , we integrate the right-hand side by parts to 
bring out a factor of e; that is, we integrate the rapidly oscillating exponential 
e lw fe*/ E and leave everything else. For the force term, we have 



f e» T (W £ , f £ ) dr = |M| £ /* eT*!^ dr 
Jo . Jo 



1 



(4.12) =e|M|^ _[^/«(t)e^+Kt/ e _^ (0 )/4(o)] 

e iw fc r / £ d , 



Here the prime on ^ indicates that terms for which u>f, = are omitted since then 
wt = 0. Introducing the integration operator 1^ defined by 



(4.14) Jo 



(4.13) LW £ (x, t) := V — <(t)Xfe- i ^ t / £ e lfe - a! , 
which is well-defined since > 1, we can write this as 

V(T, / £ ) dr = £e" ! (ir((), / £ (t)) - e (L^ £ (0), / £ (0)) 

- e f e VT {u(\ u W e , f) + (Ld* T W £ , f) + (LW £ ,d T f)} dr. 
Jo 

Similarly, integrating the next term by parts we find 

/ e UT (W e ,B{W°,W°))dT 
Jo 

, , =ee" t (LW £ ,B(W°,W o ))(t)-s(LW £ ,B(W°,W Q ))(0) 

( 4 - 15 ) r r 

e" T {v (LW £ , B(W°, W )) + %d*W e , B{W°, W )) 

+ (LW £ ,B(d T W°,W )) + (LW £ ,B(W°,d T W ))} dr. 
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Next, we consider 
ft 



I e VT (W°,B(W £ ,W £ )) dr 
Jo 

J ° jkl 

rs orsO i r>srO 

(4.16) = Ey' ^ihLZ^Ml 

jkl 3 « 



<(oW(oK(o)] 



„f rs nrsfl i nsrO 

£1 
2 



r t rs nrsO , osrO , 

/ y ^jkl+^kjl e 4( B JKW4[ m yAl dr. 



Here the prime on ^ indicates that exactly resonant terms, for which uj^ + w£ = 
and BjM+BfL® = 0, are excluded. Using the bilinear operator B^, defined for any 
W £ , W £ and W° by 



rs Di'sO i nsrO 

(4.17) (TU°,^(TUMy £ )):=^ M M ,r t& .^ e -i(o, ?+ ^)t/ 8> 

we can write (14. 16)) in the more compact form 

/ e UT (W°,B(W £ ,W £ ))dT 
Jo 

= ee» t (W°, B U {W £ , W £ )){t) - e (W°, B„(W e , W £ ))(Q) 
( 4 - 18 ) f t 

-£ / e UT {^(W ,B UJ (W £ ,W £ )) + (d T W°,B L ,(W £ ,W £ )) 



+ (w°,d;B UJ (w £ ,w e ))}dT. 

Putting these together, (|4.10p integrates to 

e ut \W E (t)\ 2 - \W E (Q)\ 2 +n f e UT \WW e \ 2 dT 

Jo 

< 2e e» l f e ) (t) -2e(LW £ , f) (0) 

(4.19) - 2e e ut (\ u W e , B(W°, W )) (t) + 2e (\ U W £ , B(W° 7 W )) (0) 

+ 2e e vt (W°, B U {W £ , W £ )) (t) - 2e (W°, B U {W £ , W £ )) (0) 



+ 2e / e^{J (r)-/i(T)+/ 2 (r)}dr. 
Jo 

Here the integrands are 

(4.20) /o := HLw £ , n + (lw £ , d T f £ ) + (\ u d;w £ j £ ), 

h := v (\ u W e ,B(W°, W°)) + (\ u d;W s , B(W°, W )) 
' 4J ' ' + {LW £ , B(d T W°, W )) + {LW £ , B(W°, d T W )), 
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(4.23) 



and 

I 2 >.= v(W 1i ,B u {y?*,W e )) + (d T W° ,B u (W e ,W E )) 

( 4 - 22 ) 

+ {W°,d* T B UJ (W £ ,W e )). 
We now bound the right-hand side of (|4.19[) . On the second line, we have 

■ *(LW e (t),r(t))-Q u W s (p),f e (p))\ 

<e^|W^(t)||/ e (t)| + |W c (0)||/ e (0)|, 
where we have used the fact that, thanks to (|3.25p , 

(4.24) \W a LW £ \<\W a W E \, for a = 0,1,2,- ■■ . 
To bound the next line, we use the estimate 

, \(W,B(W°,W))\ <C\W\ L e\W°\ L 3\VW\ L 2 

(4.25) _ , 

< C \VW\ \W°\ 1/2 \VW°\ 1/2 \VW\ 

(note that the first argument of B is W°) to obtain 

\e vt (\ u W e , B(W°, W°))(t) - %W E , B(W°, W°){0)\ 

(4.26) < e ut \VW e (t)\ \W (t)\^ 2 \VW°(t)\ 3/2 

+ |W E (0)| |W (0)| 1/2 |W°(0)| 3/2 . 

In (|4.25|) and in the rest of this proof, C and c denote generic constants which may 
not be the same each time the symbol is used; such constants may depend on M 
but not on any other parameter. Numbered constants may also depend on fi. 

We now derive a bound involving B u . Since Bjf^ + Z?jy? — in the case of exact 
resonance, we assume that uj'j + ujf, ^ 0. Then (|3.51|) implies 

(4.27) 1 m + yl <c\j\\k\. 

With this, we have for any W e , W £ and W°, 



\(W°,BUW s ,W £ ))\<lj2 



2 
jki 



r>rsO j_ TDsrQ 



KIKIK°I 



(4281 <c Y, lillfclKIKII^°l 

j+k=l 

< [ Q(x)£{x)((x) dx 3 

JJA 

<c\ww e \ L ^w e \^\w°\ Lm , 

with 1/p + 1/q + l/m = 1 and where on the penultimate line 

(4.29) 9{x) :=£;b1K|e y -",£(s) :=££|fc| K|e*'" and := e"-. 

Using (14.28)) with p = q = 2 and m — oo, plus the embedding H 2 CC L°°, we 
have the bound 

^ vt (W\B u (W £ ,W e )){t)~(W\B u (W £ ,W e ))(0)\ 



(4.30) 



< 



C (e ut \VW e {t)\ 2 \V 2 W Q {t)\ + \VW £ {0)\ 2 \V 2 W°{0)\) 
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To bound the integrand in ()4.19|) . we need estimates on dtW° and d^W £ in 
addition to those already obtained. Using the bound 

(4.31) \B°{W, W)\ L 2 < C \VW\i* < C \VW\ 2 H3/i < C \W 2 W\ 3/2 \WW\ 1/2 , 

we find from (|3.9p 

\d t W°\ L2 < C \VW\ 2 H3/i + M |V 2 ^°| + |/| 
l4 " >_) < C\\7 2 W\ 3 / 2 \\7W\^ 2 + fi\\7 2 W°\ + \f\. 

Similarly, we find from (13.341) 

\d;W s \ L2 < C \VW\ 2 H3/i + ,x \V 2 W £ \ + \f\ 
' 4 " > >! < ClV^l^lViyl 1 / 2 +n\X7 2 W £ \ + |/|. 



Now using the bound 



(4.34) \VB{W,W)\ L2 <C\V 2 W\ L1V ,\VW\ L12 < C\V 2 W\ 1 



we find 

\Vd t W°\ L2 < C \V 2 W\ 2 Hl/i + f i \V 3 W°\ + |V/°| 

< C \V 3 W\ 1/2 \V 2 W\ 3/2 + /x |v 3 iy°| + |V/°|, 

(4.35) 

\sjd* t w £ \ L2 < c \x/ 2 w\ 2 H1/4 + n \v 3 w £ \ + |v/ e | 

< C |V 3 W| 1/2 |V 2 W| 3/2 + fi \W 3 W £ \ + |V/ £ |. 
The bound for Jo follows by using (|4.33l) . 

\i \ L , < c (\vw\ 2 H3/i + (n + c) \v 2 w £ \ + \r\) (in + |9 t r i) 

< C (\VW\ 2 H3/i + (ju + c) |V 2 JU| + ||/|U ||/|| g , 

where we have used the fact that \V a W £ \ 2 < \V a W £ \ 2 + |V Q iy°| 2 = |V Q Fy| 2 . 
Next, using (|4.28|) we bound 1% as 

\h\ L 2 < fic\W°\ Lao \WW £ \ 2 + c\d T W°\ \WW £ \ 2 L , 

+ c\W°\ Loa \VW £ \ \Vd*W £ \ 

< ^ |V 2 M/°| |VT^-| 2 + c (\VW\% 3/4 + n \V 2 W°\ + |/°|)|V^| 2 ff3/1 
,4 " > ' 1 +c|V 2 iy°||Viy £ | (\V 2 W\ 2 H1/i + M |V 3 V^ £ | + |V/ £ |) 

<c\VW\ \V 2 W\ |V 2 IT^ 1/4 + fic \V 3 W\ \V 2 W\ |VW| 
+ |V 2 V7| 3/2 |VTU| 1/2 ||/|| g 

where interpolation inequalities have been used for the last step. The bound for I\ 
is majorised by that for 12- 

Putting everything together, we have from (|4.19p 

e vt \W £ (t)\ 2 - \W £ {Q)\ 2 

< ec 2 e vt |V 2 W^)|(|W(t)| 2 + ||/y + £C2 |V 2 Wo|(|VWo| 2 + ||/|| g ) 

(4.38) r* 

+ ec 3 / e^{\W\ m \W\ H2 \W\ 2 Ha/i + ^\W\ H3 \W\ H2 \W\ H1 
Jo 

+ (\W\ 3 ^\W\^ + (/i + c) \W\ H2 + ||/|| g ) \\f\\ g } dr. 
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Now by flUEl) and ([219]) . we can find K*(||/|| g ) and T*(|VW |, ||/|| g ) such that, 
for^>T», 

(4.39) c |V s ^(t)| 2 + (ji + c') \W s W(t)\ + ||/|| g < K* 

for s £ {0, 1, 2, 3}. Let t' :=t- T„ and relabel t in (|Q5|l as t'. We can then bound 
the integral in (|4.38p as 

(4.40) f e" T {---}dr< '—— ± c 4 K*(\\f\\ g ) 2 . 
Jo v 

Bounding the remaining terms in (|4.38|) similarly, we find 
(XA ' ' < e -"(*- T *) |^(T*)| 2 +£c 5 (X 2 + ^ 3/2 ). 

This proves the theorem, with # g (||/|| g ) 2 = 2 c 5 (KI+kI /2 ) and T g (|VW |, ||/|| g , e) = 
T.-logfecs (X,+Xy 2 )]/^. 

5. Higher- Order Estimates 

When dtf = in the very simple model (|4.3p . we can obtain a better estimate 
on i' = x — U where U — ef /(e/i + i) than on x, namely that x'(t) — > as 
t — * oo; here U is the (exact, higher-order) slow manifold. The situation is more 
complicated when / is time-dependent, or when x is coupled to a slow variable y 
with the evolution equations having nonlinear terms. In this case, it is not generally 
possible to find U (explicit examples are known where no such U exists), and thus 
x'(t) as t — > oo for any U(y,f;e). Nevertheless, it is often possible to find a 
U* that gives an exponentially small bound on x'(t) for large t. We shall do this 
for the primitive equations. 

More concretely, in this section we show that, with reasonable regularity as- 
sumptions on the forcing /, the leading-order estimate on the fast variable W e 
in the previous section can be sharpened to an exponential-order estimate on 
W e - U*{W°J;s), where U* is computed below. As in [31], we make use of 
the Gevrey regularity of the solution and work with a finite-dimensional truncation 
of the system, whose description now follows. 

Given a fixed k > 0, we define the low-mode truncation of W by 

a 

(5.1) W<(x,t) = (P < W)(x,t) := ti£Xjg , e- iw 2*/ e e i *-" 

\k\<K 

where the sum is taken over a 6 {0, ±1} and k £ with |fc| < n. We also define 
the high-mode part of W by W > :— W — W < . The low- and high-mode parts of 
the slow and fast variables, W 0< , W 0> , W e< and W e> , are defined in the obvious 
manner, i.e. W a< with a — in (|5.1[) and W e< with a £ {±1}. It is clear from 
(|5.1[) and (|3.18[) that the projection P < is orthogonal in H s , so P K commutes with 
both A and L in (f2"3T]) . We denote P < BbyB < . 

It follows from the definition that the low-mode part W K satisfies a "reverse 
Poincarc" inequality, i.e. for any s > 0, 

(5.2) \VW<\ H . <k\W<\ h .. 
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If W G G" T (M), the exponential decay of its Fourier coefficients implies that W > is 
exponentially small, that is, for any s > 0, 

(5.3) \W>\ H . <CU s e-™|^| G „. 

The first inequality evidently also applies to the slow and fast parts separately, i.e. 
with W K replaced by W 0< or W £< ; as for (|5.3|) . it also holds when W > on the lhs 
is replaced by W 0> or W £> . 

We recall that the global regularity results of Theorem imply that, with Gevrey 
forcing, any solution W G H 1 (l\i) will be in G" T (M) after a short time. As in [31] 
and following [22], the central idea here is to split W e into its low- and high- mode 
parts. The high-mode part W e> is exponentially small by (|5.3p . We then compute 
U*(W 0< , / < ; e) such that W £< — U* becomes exponentially small after some time. 

Following historical precedent in the geophysical literature, it is natural to present 
our results in two parts, first locally in time and second globally. (Here "local in 
time" is used in a sense similar to "local truncation error" in numerical analy- 
sis, giving a bound on the time derivative of some "error".) The following lemma 
states that, in a suitable finite-dimensional space, we can find a "slow manifold" 
W e< = U*(W 0< , / < ; e) on which the normal velocity of W e< is at most exponen- 
tially small: 

Lemma 2. Let s > 3/2 and rj > be fixed. Given W° E H S (M) and f € H S (M) 
with dtf = 0, there exists e**(\W Ifln^tV) such that for e < £*» one can find 
k(s) and U*(W , f < ;e) that makes the remainder function 

n*(W Q< J < ;e):=P < [(DU*)G*} + -LU* 

(5.4) e 

+ B £< (W 0< + U*, W 0< + U*) + AU* - f e< 

exponentially small in e, 

(5.5) \n*(W°<J<;e)\ H , <c r [(\W°<\ HS + \f\ HS ] expHy/e 1 / 4 ); 
here DU* is the derivative of U* with respect to W 0< and 

(5.6) G* := -B Q< (W 0< + U*,W 0< + U*) - AW 0< + f a< . 

Remarks. 

1. The bounds may depend on s, /i and M as well as on rj, but only the latter is 
indicated explicitly here and in the proof below. 

2. Given k fixed, U* lives in the same space as W £< , that is, (W°, U*) L2 = and 
P < U* = U*. 

3. In the leading-order case of <g] the slow manifold is J7° = and the local error 
estimate is incorporated directly into the proof of Theorem [TJ we therefore did not 
put these into a separate lemma. 

4. Unlike formal constructions in the geophysical literature (see, e.g., [5j [37]), 
our slow manifold is not defined for all possible W° and e. Instead, given that 
IW-^Ig^ < R, we can define U* for all e < a); generally, the larger the set of 
W° over which U* is to be defined, the smaller e will have to be. 

5. In what follows, we will often write U*(W°J;e) for U*(P < W°, P < /; e); this 
should not cause any confusion. 
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Using the Lemma and a technique similar to that used to prove Theorem [TJ we 
can bound the "net forcing" on W — W £< — U* by 1Z*. The dissipation term AW 
then ensures that W' eventually decays to an exponentially small size. This gives 
us our global result: 

Theorem 2. Let W £ if^M) and V/ € G" T (M) be given with d t f = 0. Then 
there exist £*(/;<r) and T*(|VWo|, |V/| G „) such that for e < e* and t > T*, we 
can approximate the fast variable W £ (t) by a function U*(W°(t), /; e) of the slow 
variable W°(t) up to an exponential accuracy, 

(5.7) \W s (t)-U*(W°(t),f;e)\ L2 <K*(\Vf\ G <r,o-) exp^a/e 1 / 4 ). 

As in Theorem [1] here is a continuous increasing function of its arguments; 
W(t) = W°(t) + W e {t) is the solution of (pTTTTl with initial condition W(0) = W . 
As before, the bounds depend on /i and M, but these are not indicated explicitly. 

Remarks. 

6. With very minor changes in the proof of Theorem [2] below, one could also show 
that, if / e H n+1 and d t f = 0, then \W £ (t) - U n (W°(t), f;e)\ L2 is bounded as 
e ™/ 4 f or sufficiently large n and possibly something better for smaller n. 

7. Recalling remark [¥] above, our slow manifold is only defined for e sufficiently 
small for a given |VK 0< | (or equivalently, for |W /0< | sufficiently small for a given e). 
The results of Theorem tell us that Wit) will be inside a ball in G" T (M) after 
a sufficiently large t; we use (twice) the radius of this absorbing ball to fix the 
restriction on e. Thus our approach sheds no light on the analogous problem in the 
inviscid case, which has no absorbing set. 

8. As proved in [121 E2 H3] , assuming sufficiently smooth forcing, the primitive 
equations admit a finite-dimensional global attactor. Theorem [2] states that, for 
£ < £*(|/Ig")i solution will enter, and remain in, an exponentially thin neigh- 
bourhood of U*(W Q< , / < ;e) in L 2 (M) after some time. It follows that the global 
attractor must then be contained in this exponentially thin neighbourhood as well. 

9. The dynamics on this attractor is generally thought to be chaotic [30]. Thus 
our present results do not qualitatively affect the finite-time predictability estimate 
of[S]. 

10. When dtf ^ 0, the slaving relation U* would have a non-local dependence 
on t. Quasi-periodic forcing, however, can be handled by introducing an auxiliary 
variable 6 = (9i, ■ ■ ■ ,0 n ), where n is the number of independent frequencies of /. 
The slaving relation U* would then depend on as well as on W a< . 

11. Bounds of this type are only available for the fast variable W £ ; no special 
bounds exist for the slow variable W° except in special cases, such as when the 
forcing / is completely fast, (W°,f) L2 = 0. 

We next present the proofs of Lemma [5] and Theorem O The first one follows 
closely that in [31] which used a slightly different notation; we redo it here for 
notational coherence and since some estimates in it are needed in the proof of 
Theorem [2] As before, we write (•, •) = (•, -)^2 and | • | = | • \ L2 when there is no 
ambiguity. 
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5.1. Proof of Lemma [2} As usual, we use c to denote a generic constant which 
may not be the same each time it appears. Constants may depend on s and the 
domain M (and also on fi for non-generic ones), but dependence on r\ is indicated 
explicitly. Since s > 3/2, H S (M) is a Banach algebra, so if u and v G H s , 

(5.8) \uv\ s < c\u\ s \v\ s 

where here and henceforth | • \ s := | • \ Hs . Let us take s < 1 and k as given for now; 
restrictions on e will be stated as we go along and k will be fixed in (15.22|) below. 
We construct the function U* iteratively as follows. First, let 

(5.9) -LU 1 = ~B £< (W a< : W Q< ) + f £< , 

e 

where U 1 € range L for uniqueness; similarly, U n S range L in what follows. For 
n = 1, 2, • ■ • , let 

(5.10) -LU n+1 = -P < \(DU n )g n ] ~B £< (W 0< + U n ,W 0< +U n )-AU n + f £< , 
e 

where DC/™ is the Frechet derivative of U n with respect to W 0< (regarded as living 
in an appropriate Hilbert space) and 

(5.11) g n := -B 0< (W a< + U n , W 0< + £7") - AW 0< + f<. 

We note that the right-hand sides of (|5.9p and (|5.10p do not lie in kerL, so U 1 
and U n+1 are well defined. Moreover, U n lives in the same space as W £< , that is, 
U n e P < range L; in other words, (W°, U n ) = and P < U n = U n . 

For 7] > 0, let D ri (W 0< ) be the complex ry-neig hbourhood of W 0< in P iP(M). 
With W 0< defined by (EH), this is 

D V (W°) = lw° : W°(x,t) = ™l X k eikX with 

|*l<« 

(5.12) 

_ fc3) and J2 |*| Sta |ta£-t«£| a <ir , |. 

|fc|<« * 

Since W°(x,t) and X% are real, u>£ must satisfy (|3.29a ). but in (|5.12[) need 
not satisfy this condition although it must satisfy (|3.29b V We can thus regard 
D V (W 0< ) C {(w k ) : < |fe| < k and w {klM k ) = u> fclifc2:fc3 )} = C m for some m. 
Let 5 > be given; it will be fixed below in (|5.22p . For any function g of W a< , let 

(5.13) \g(W°<)\ s . n := sup | 5 (W)| S ; 

this expression is meaningful when D v - n $(W 0< ) is non-empty, that is, for n € 
{0, • • • , L^/^J =: For future reference, we note that 

(5.14) \W a <\ s . fi <\W°<\ s +r h 

Our first step is to obtain by induction a couple of uniform bounds (|5.25p - (|5.26[) . 
valid for n € {1, • • • ,«*}, which will be useful later. First, for U 1 , we have 

(5.15) \\LUX, < \B £ <(W°<,W°<)\ s . A + \f £ <\ s 

which, using the estimate \B(W, W)\ s < c\VW\ 2 s and (|5.2[) . implies 

(5.16) lUX.Ksc^K'lW^ + m). 



W (ki,k 2 ,k 3 ) ~ W (kiM 
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Next, we derive an iterative estimate for |J7™| S .„. Using the fact that |-| s . m < 
whenever m > n, we have for n — 1, 2, • • • , 

^ 7 |C/" +1 L;n + l < \(OU n )G n \ S[n+1 + \B S <(W°< + IT\ W°< + U n )\ S ; n 

(5.17) £ 

+^ 2 \w°% n + \r<\ s . 

The first term on the right-hand side can be bounded by a technique based on 
Cauchy's integral formula: Let D v (z ) C C be the complex 77-neighbourhood of z . 
For ip : D n (z ) — > C analytic and 6 £ (0, 77) , we can bound \ip'\ in D^s^Zq) by \tp\ 
in D v (z ) as 

(5.18) y ■ z\ Dri _ s{Zo) < g\(p\ Dr>(Zo) \z\ c . 

Now by (|5.9[) t/ 1 is an analytic function of the finite-dimensional variable W 0< , so 
assuming that U n is analytic in W 0< we can regard the Frechet derivative DU n as 
an ordinary derivative. Taking for (p' in ()5.18j) the derivative of U n in the direction 
Q n (i.e. working on the complex plane containing and Q n ), we have 

(5.19) \(DU n )g n \ s . n+ l<^\U n Un\g n \ s -,n- 

Using the estimate 

(5.20) \B £< {W 0< + U n ,W 0< +U n )\ s . n < c\V(W 0< +U n )\% n < ck 2 \W 0< + U n \% n 
we have 

\U n+ % n+1 < j \U n \ s , n (cn 2 \W°< + U% n + fiK 2 \W°<\ s . n + \f<\ s ) 
+ en 2 c \W 0< + U n \l. n + fien 2 \U n \ s - n + e |/ £< | s ■ 
To complete the inductive step, let us now set 
(5.22) S = e 1/4 and k = £~ 1/4 . 
With this, we have from (|5.21[) 

\u n+1 \ s , n+1 < £ 1/4 C! \u n \ s .J\w°< + u% n + v\w°<\ s:n + e v 2 |/ 0< L) 

+ s 1/2 <h (|VF°< + UXn + V W n Un + e 1/2 \f<\ s ). 
We require e to be such that 

(5.24) e 1 / 4 (c + c 1+ c 2 ) {\W°% + fi \W°<\ S:0 + \f\ s ) < \ min{l, \W°<\ S } 
and claim that with this we have 

(5.25) \U n \ S ;n < ^ 1/4 % (\W°<\ 2 S , Q +/i|VF°<| s;0 + |/<| s ) 

with c v — 4 (c + c x + c 2 ). Now since e < 1, (|5.16|) implies that it holds for n = 1, 
so let us suppose that it holds for m = 0, ■ • • , n for some n < n*. Now (15.24)1 and 
([5T23| imply that 

(5.26) |£H s;m <|TU <L<|TU°<U and \U m \ s;m < 1 
for m = 0, • • • Using these in (|5.23[> . we have 

|c/" +1 | s; „ + i<4 e 1 / 4 Cl (|w r0< l^ ; o + M^ 0< L ;0 + |/ < L)|c/"| s; „ 

(5-27) + ^ 2 c 2 {\W <\% + f x\W <\ s . fi + \f<\ s ) 

< (|T^°<lf ;0 + / x|W°<| s;0 + |/<U). 
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This proves (|5T25]) and (IQ6]) for n = 0, • • • ,n*. 
We now turn to the remainder 

(5.28) U° :=B e< (W 0< ,W 0< ) - f s< 
and, for n = 1, • • • , 

(5.29) TZ n := P < \(DU n ) G n ] + -LU n + B e< {W 0< + U n , W 0< + U n ) + AU n - f e< . 

e 

We seek to show that, for n = 0, • • • ,n*, it scales as e _Tl . We first note that by 
construction TZ n ^ ker L, so L _1 7?. n is well-defined. Taking U° — 0, we have 

(5.30) n n = jL(U n -u n+1 ). 

We then compute 

n n+i = p < [(DC7"+ 1 )^«+ 1 ] + IlC/" +1 

iS<nx/0< i TT n + 1 M/0< i ZT™+1\ _i_ ATTn+l 



(5.31) 



+ B e< (W u< + U n+i ,W"< + U n+L ) + AU n+i - f s< 

= p < [(Du n+1 )(g n + sg n )] + -lu 71 - n n 

e 

+ B e< (W 0< + U n , W 0< + U n )-e B e< (W 0< + U n , L~ lr R. n ) 



~eB £< {L- l n n , W 0< + U n+1 ) + AU n -eAlr x Tl n - f E< 
= P < [(DU n )Sg n ]-eL- 1 P < [(DTZ n )g n+1 ]-eAL- 1 TZ n 

-eB e< {L- lr R n , W 0< + U n+l ) -eB £< (W 0< + U , \L- 1 1l n ) 1 
where we have used (|5.30|) and where 

5 gn ._ gn+1 _ gn 

(5 32) 

= sB 0< {W 0< + U n+ \L- 1 K n )+eB 0< (L- 1 K n ,W 0< +U n ). 
To obtain a bound on JZ n , we compute using (|5.26|) 

\g n \ s;n < c (MW°< + U n )\%n + M |A^ 0< L ; „ + |/ 0< L) 

( ■ j <c K mw°<\i fi +v\w°% 0+ \f\ s ), 

as well as 

|^"L, n+1 < ec|V(^°< + C/"+ 1 )| s; „ +1 |VL- 1 ^"| s; „ +1 

(5.34) + ec |VL- 1 ft n | Sin+1 |V(W°< + CP)|,.„ 

<e K 2 c\^\ n+1 \W°<\ s . fi . 

(Note that we can only estimate Sg n in D , n _^^ s and not in D v _ nS ; similarly, 
since the definition of 1Z n involves DC/", it can only be estimated in D , n+1 ^ g .) 
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Wc then have 

l^ n+ N s; „ +2 <|D^ n | s;n+1 |^"| s; „ +1 + £ |i- 1 D7e"| s . n+2 |g™+ 1 | s;n+1 

+ \K% n+1 + e |VL- 1 ft n | s;n+1 |V(W < + U n+1 )\ S]n+1 
+ e\V(W°<+U n )\ s .jVL~ 1 K n \ S ; n+1 



(5.35) i 



< -\U n \ s , n e n 2 |^"| s; „ +1 |^°<| s;0 + c | |^"| s;n+1 |C?" +1 | s; „ +1 

+ Ask 2 \K n \ s;n+1 \W°<\ s;0 + SfiK 2 \1Z n \ s;n+l 

< e 1 /* n-, n+ i c e (\w 0< \l + n \w°<\ s . fi + |/<| s + M ) 

where for the last inequality we have assumed that 

(5.36) e 1 / 4 < min{ /i /|W 0< | s . ,M%/4}. 

If we require e to satisfy, in addition to e < 1, (|5.24[) and (|5.36[) . 

(5.37) e 1 / i c e (\W <\ 2 . + f i\W % + \f<\ s + M ) < ± , 
we have, for n = 0, 1, • • • , n* — 1, 

(5.38) |^ +1 L. ;n+2 <i|^| s;n+1 . 
Along with the estimate 

(5.39) IftV^flW^fto + im 
taking rt = rt* — 1 leads us to 

, c . I^V < I^^U, < c r (|W 0< |f ;0 + |/<|.) cM-ri/s 1/4 ) 
(5.40) 

< Cl ,[(|M/°<L + ?7 ) 2 + |/<L] eM-v/e 1/4 )- 

The lemma follows by setting U* — U n *~ 1 and taking as e** the largest value that 
satisfies e < 1, (pT2Ij) . (jOBl) and (|Q7|l . 

For use later in the proof of Theorem [21 we also bound 

|V(1 - ? < )[(DU*)g*]\ L2 < ce-™\{w*)g*\ 2 , nt 
( 5 - 41 ) <ce—\\u*\ 2 ^ 1 \g*\ 2 ^ l 

<ce-™K 2 (\W°<\l + v\W°<\ 2 ; + \f\ 2 ) 2 

where for the last inequality we have used (|5.25| and (|5.33|) with n = n* — 1. 

5.2. Proof of Theorem [2j We follow the conventions of the proofs of Theorem [T] 
and Lemma [2] on constants. We will be rather terse in parts of this proof which 
mirror a development in the proof of Theorem [1] 

First, we recall Theorem and consider t>T:~ max{T 2 , T a } so that | V 2 W(t) | < 
K 2 and |V 2 V7(i)| G „ < M a . We use Lemma[5]with s = 2 and, collecting the con- 
straints on e there, require that 

e 1 / 4 c v ((K 2 + V ) 2 + /z (K 2 + v ) + |/| 2 ) < | min{l, K 2 }, 
(5.42) £l,A < min {M/(^2 + ?7),MC(7/4, 1}, 

e 1 '* c e ((K 2 + 7 1 ) 2 +^(K 2 +r 1 )+^+ |/| 2 ) < 
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where c e is that in (15.37p . (We note that all these constraints are convex in K2, so 
they do not cause problems when |M /0< | < K%.) Further constraints on e will be 
imposed below. We note the bound (|5.40[) and 

(5.43) \U*\ H2 < \U*\ 2 . Ht < £ 1/4 c v {(K 2 + r\) 2 + (i (K 2 + v ) + |/| 2 ) 

which follows from (|5.26p . 

We fix k = e -1 / 4 as in (|5.22[) and consider the equation of motion for the low 
modes W < , 

d t W < + -LW< + B < (W < , W<) + AW< - f< 
e 

(5.44) = -B < {W > ,W)- B < {W < ,W > ) 

=:H. 

Writing 

(5.45) W £< = U*{W 0< J<;e) + W , 

the equation governing the finite-dimensional variable W'(t) is 

d t W + -LW' + B £< {W < ,W < ) + AW' 

(5.46) £ 

= -d t U* - -LU* - AU* + f e< + W. 

£ 

Using (|5.4p . this can be written as 

d t W + -LW + B e< {W < , W) + B £< (W', W 0< + U*) + AW' 

e 

(5-47) = -K* - {l-P < )[{DU*)g*}+H £ 

=: -TV +H e . 
Multiplying by W' in L 2 (M), we find 

(5.48) ~\W'\ 2 + (W',B £< (W',W 0< +U*))+fi\VW'\ 2 = ~{W' ,K*) + [W' ,H e ). 

We now write the nonlinear term as 

(W',B e< (W',W 0< + U*)) = {W',B(W',W 0< + U*)) 

(5.49) = {W, B(W', U*)) + [W', B(W', W 0< )) 

= (W, B(W', U*)) - {W 0< , B{W', W')). 

Following the proof of Theorem [1] [cf. (|4.10[) ]. we rewrite (|5.48|) as 

(5 50) ^(^'l^'l 2 ) +M e,y *|VPF'| 2 < -2e vt (W \TV) +2e ut [W \W) 

- 2 e vt (W, B(W', U*)) + 2 e vt (W 0< , B(W', W')). 

We bound the first two terms on the right-hand side as 

2\{W',H*)\ < -\VW'\ 2 + -\K*\ 2 , 
(5.51) 6 M 

2\(W',H £ )\<-\VW'\ 2 + -\H £ \ 2 . 
6 fi 
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(5.56) 



As for the third term in (|5.50[> . we bound it as 

2 \(W, B(W, U*))\ < c \W'\ L e\VW'\ L 2 \VU*\ L i 

' "' ; ' 2 ' < I Wf Cl £ V4 (( | W 0< |a + v) 2 + ^ + v) + |/<|a) 

where we have used (I5.43[) in the last step. We now require e to be small enough 
so that 

(5.53) e 1 / 4 c 1 ((K 2 +f 1 ) 2 + ^K 2 + f i V +\f\ 2 ) < |, 
which implies that, since |M /0< | 2 < by hypothesis, 

(5.54) 2\{W',B(W',U*))\ < -\VW'\ 2 . 

6 

With these estimates, (|5.50[) becomes 

(5 55) d7 (e " |W " |2) + 2 6 " |Wf " £ 6 " (l ^' 2 + |H£|2) 

+ 2e vt (W 0< ,B(W',W')). 

Integrating this inequality and multiplying by e _l/T , we find 

11 r T+t 

e ut \W'(T + t)\ 2 -\W'(T)\ 2 + ^ / e^-^lWfdr 

2 7t 

< r +f e" (T - T) |-(|^T + |H e | 2 ) + 2(^ 0< ,B(W r ',iy'))} dr. 
We then integrate the last term by parts as in (|4.18[) . 

y + e" (r - T) (^ 0< ,S(W^',T^'))dT 

= ee 1 ' 4 (VK 0< , B U (W, W')){T + 1) - e (W 0< , B U (W', W')){T) 

(5-57) 

-e / e^ T - T ){^(^ 0< ,^(^',Ty')) + (a T iy 0< ,^(^',^')) 

+ 2(vK 0< ,B w (a;v^',^'))}dT. 

To bound the terms in the integral, we first need to estimate 

\VB £< (W',W a< + U*)\ L2 < k\B s< (W',W 0< + U*)\ l2 

<ck\VW'\ l2 \V(W°<+U*)\l~ 
<ck 2 \VW'\\W 0< +U*\ H 2 

<ck 2 \VW'\\V 2 W°\ 

where for the last inequality we have used (|5.26p . Using this and the bound 

(5.59) \VB e< (W < ,W')\ L3 < cK|VW r< | i0 o|Viy'| L 2 < ck 2 \V 2 W\ |VW'| 
for the term B £< (W<, W) in (|CT7|) gives us 

(5.60) \VdlW'\ L 2 < ck 2 \VW'\ \V 2 W\ +fiK 2 \VW'\ + \VK*\ + \VH £ \. 



26 TEMAM AND WIROSOETISNO 

The worst term in (|5.57j) can now be bounded as 

e\(W 0< ,BUd;W',W'))\ <ec\W 0< \ Loo \VW'\ L2 \Vd;w'\ L2 
(5 61) < c 2 e k 2 K 2 |VW"| 2 + c 3 e n 2 M K 2 \VW'\ 2 

+ ^ |W| 2 + (\vn*\ 2 + |vh £ | 2 ). 

If we now require that e satisfy 

(5.62) e x ^c^Kl<^ and £ V* C3 K 2 <±, 

we have 

cirw 0< _ n..(a*w'.w'w\ < 

16 1 1 ii 



(5.63) e|(w°< s^w.W))! < -^\yw'\ 2 + ^^-(\vn*\ 2 + \vn £ 



(5.66) T 



Bounding another term in (I5.57[) as 

e\(dtW Q< ,B^W',W'))\<ec\dtW°\ L ^\VW'\ 2 L2 

(5.64) < ec\d t W°\ H2 \WW'\ 2 L2 

< ec 4 ( K K 2 + fiK 2 K 2 + \ f\ 2 )\WW'\ 2 

and requiring that e also satisfy 

(5.65) eVa^^ + ^ + i/i^iL 

plus a similar estimate for the first (easiest) term in (|5.57[) . we can bound the 
integral on the r.h.s. as 

pT+t 

/ \(W 0< ,B(W',W))\dT 

<H [ + V( r ~ T )|Wf dT + ^if 2 (||V^ 2 + ||Wi e || 2 )(e" 4 -l) 
2 Jt M 

where ||Vft*|| := sup^o,^ \\7K*{W a , f;e)\ and similarly for ||VW S ||. Bounding 

the limit term in (|5 . 57[) as 

(5.67) \(W Q< ,B^W',W'))\<C\W Q< \ L ^\VW'\ 2 L2 <cK 2K 2 \W'\ 2 , 
(|5.56p becomes 

(l-e^ 2 c 5 K 2 )\W'(T + t)\ 2 

(5 ' 68) < e~ vt (1 + e l ' 2 c b K 2 ) \W\T)\ 2 + ^ K 2 2 {\\VR*\\ 2 + \\VW\\ 2 ). 

M 

To estimate ||VH e ||, we use (jiQ4"l) . (|53|) and (|2~20|) . to obtain 

«™ |V5 e< (W>, W)\ L2 + \X7B £ <(W<, W>)\ L , < CK\VW>\ L <\VW\ Lt 
(5.69 

<cne-' yK M a K 2 . 

Now (|5.41|) implies that 

(5.70) | V(l - P < )[(DJ7*)£*]| L2 < ce"™ k 2 ((tf a + V f + M (Jf 2 + „) + |/| 2 ) 2 ; 
this and the previous estimate give us 

(5.71) \\VH e \\ L 2 < ce- k 2 [M a K 2 + {{K 2 + ry) 2 + fx (K 2 + r?) + \f\ 2 ) 2 ] . 
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Meanwhile, using (|5.40[) we have 

(5.72) \\VR*\\ L2 < c ((K 2 + nf + |/| 2 ) cxp(- V /e^). 

Setting 7] = o and requiring e to satisfy, in addition to (|5.42|) , (|5.62|) and ()5.65|) , 



(5.73) £ 1/2 c 5 if 2 <i, 
we have 



\W\T + t)\ 2 <4e- ut \W'(T)\ 2 
(574) + j 2 [(K 2 + of + ^(K 2 + of + l/l! + + M 2 ^ 2 ] exp(-2a/ £ 1 /4 ) . 

Since |W'(T)| < ci^2 by Theorem 0, by taking t sufficiently large we have 

\W'(T + t)\<- [{K 2 + of + M 2 (X 2 + of + |/| 2 + \f\i + M 2 a K 2 2 \ x 

(5.75) M 

e 1 / 2 exp(-a/e 1 / 4 ). 

And since 

(5.76) \W e -U*\ 2 < \W e> \ 2 + \W'\ 2 <cM 2 exp(-2a/e 1 / 4 ) + |^'| 2 , 
the theorem follows by the same argument used to obtain Theorem [TJ 

Appendix A. 

Proof of Lemma [1} Since B^ — when jzk 3 \l\ — 0, we assume that |Z| 7^ 
in the rest of this proof. As before, all wavevectors are understood to live in Z^ — {0} 
and their third component take values in {0, ±27r/Z/3, ±47r/L3, • • • }. 

We start by noting that an exact resonance is only possible when j and fc lie on 
the same "resonance cone", that is, when |j|/|j3l = 1^1/1^3 1, or equivalently, when 
Ij'l/lisl = 1^' l/l ^3 1- There are only two cases to consider: 

(a') When f = k' = 0, we have Bffi = Bf$ = 0. 

(b') In the generic case jzk^j'] \k'\ ^ 0, direct computation using the resonance 
relation r\j\/j 3 + s|fc|/fc 3 = gives Bjf^ + = 0. This result also follows as the 
special case u>j + w£ — in (|A.9|) below. 

Now we turn to near resonances. There are several cases to consider, and we 
start with the generic (and hardest) one. 

(a) Suppose that \j'\ |fe'| f with V f 0. We define Q and 9 by 

(A.l) 29. := w r 3 - w% and 20fi := + 

(We note that f2 and 9 could take either sign. Our concern is obviously with small 
\9\, when when uij and w£ are nearly resonant, so we will restrict 9 below.) Now 
this implies that 

(A.2) w r 3 = (1 + 0)fi and w%=(d- 1)0. 
We first note that 

(A.3) |j'| 2 /|j 3 | 2 = (l + 0) 2 tt 2 -l and |fe'| 2 /|fe 3 | 2 = (1 - 9) 2 Q 2 - 1 
and compute 

(A.4) |j'| 2 ^ - |fe'| 2 ^ = 49n 2 j 3 k 3 = --*L rs \j\ \k\. 

J3 k 3 1 - 9- 
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Direct computation gives us 

r>rsO i nsrO 
n jkl ' n kjl 

- IottwwI^' - «"> + (p ' + P " )Q + PV + p "«"i 

where 

P := f A k' Q := -j' ■ k 1 

, i r \j\ , , .r\j\k 3 

J3 J3 J3 

p« :=i flM (j -^ fe0 _i£iM|= 

K 3 h k; 

After some computation, we find 

p> + p" = 2en i (j> • k' + J^yi | fe | - |J "' 2 * 3 |fc '' 2 J3 



(A.6) 



bT 


Q' : 


«3 


«3 


-Ifc'l 2 

3 


Q" 


J3 


^3 | . 
J3 



1-02 I'M I 2 J3 2 fc 3 

Q' - Q" = 2^ |j| |fe| - i(j' A fe')) , 

(A.7) V 1 61 _ y 

P'Q' + P"Q" = 2 en\ 2 \ sljl W r M s M (j > A fe') + 2ir S |j| |fc|(j' • fe') 

l i - e z j 3 fc 3 

+ ^(/-fe')(|fc'l 2 | i + b-'l 2 ^)-ib'l 2 Ifc'l 2 }, 



from which we obtain 

(A 8) + fcj7 - 2jij t 1 (J ' } bib'iifei ifc'i - 21 uiifei 

2irs0 2 (j' -k')j 3 k 3 2rs/nj'Ak' 2/V \j\ \k\ , } 

1 6* \f\ |fc'| + 1 - 6* \j'\ |fc'| J3 " 3 + l-02 \j'\ |fc'| " * J j • 
Now if we require that 1 6* | < 6*o < 1 , we have the bound 

(A.9) + flgO| < M (4 + T A f ) ^ K + -||. 

To take care of the case |6*| > 9q, we note that in this case 

(a.io) | wJ+wi |>«b(|L + M). 

We note that since #o < 1 by hypothesis, this inequality holds both when u)jU>^ < 
and ojjWk > 0. Using (|3.49p , we then find 

(A.ll) \B^ l+ B^\ < \B^\ + \B^\ <V5|M|(|fe'| + |i'| + |fe'|jg| + |j'|||| 
Putting these together, we find after a short computation, 
(A.12) |B$ + BJ5»| < ^^(\js\ + |fc3|)K + 
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(b) Suppose now that \j'\ \k'\ ^ but l' = 0. Wc find using j' + k' = 0, 

(A.13) + Bj$ = i |M| 5 j+k - t ~ 1Sg 2 ^||fc| |fc/| (i3 + h*M + u' h ), 
and thus the bound 

(A.14) |^+Bj^|<^(|j 3 | + N)K + a;i|. 

(c) Finally, we consider the case j = and fc' ^ (which obviously implies the 
case k = and j 7^ 0). After some computation using Z' = fc', we find 

(A.15) 2$? + BSS? = ^^ hih ~ irk 2 )\k'\ (sr - £ 
But since in this case 

(A.16) |wj = |rsgnj 3 -s|fe|/fc 3 | = \rs - \k\/k 3 \, 
we have the bound 

which holds whether or not Z3 = 0. We recall that there is nothing to do when 
f = k' = since then = Bf$ = 0. 

The lemma follows upon fixing 6* and collecting (|A.9|) ■ (|A.12[> . (|A.14|) and (|A.17j) . 
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